Elliptic and Triangular flow in asymmetric heavy-ion collisions 
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We present a study of the elliptic (112) and triangular (173) flow and their corresponding eccentricity 
fluctuations for asymmetric (Au+Ag, Au+Cu and Au+Si) coUisions at ^/sm^ = 200 GeV. These are 
compared to the corresponding results from symmetric (Au+Au and Cu+Cu) collisions at the same 
energy. The study which is carried out using a multi-phase transport (AMPT) model shows that 
triangularity (es), fluctuations in triangularity and V3 do not show much variation for the different 
colliding ion sizes studied. However the eccentricity (£2), fluctuations in eccentricity and «2 shows 
a strong dependence on colliding ion size for a given number of participating nucleons. Our study 
thus indicates that asymmetric heavy-ion collisions could be used to constrain models dealing with 
flow fluctuations in heavy-ion collisions. 

PACS numbers: 25.75.Ld 



I. INTRODUCTION 



Knowing the initial geometry and fluctuations in 
heavy-ion collisions has recently been shown to have 
important consequences on interpretating the data 
from the Relativistic Heavy Ion Collider (RHIC)and 
the Large Hadron Collider (LHC) experiments 
For example, the long standing physical interpreta- 
tion of dihadron correlations in azimuthal angle Q 
and pseudorapidity [1) to be due to propagation of 
shock waves (due to high momentum quarks or glu- 
ons) in the medium and jet-medium interactions are 
being revisited [l|, Q . The contribution from the odd 
harmonics associated with the particle azimuthal an- 
gle distribution (originally thought to be zero due 
the left-right symmetry in the transverse plane of 
symmetric heavy-ion collisions) to dihadron correla- 
tions are found to be important. Now it is being 
widely discussed that non zero odd harmonic con- 
tributions could arise from fluctuations in the trans- 
verse positions of nucleons undergoing hadronic scat- 
tering. Further it is now known that even the ideal 
hydrodynamics calculations under predicts the mea- 
sured elliptic flow in central heavy-ion collisions at 
RHIC [6]. Only inclusion of flow fluctuations aris- 
ing due to the eccentricity fluctuations can provide 
some satisfactory explanation of the data. It is also 
believed that measurement of higher order flow coef- 
ficients and their fluctuations can substantially im- 
prove the constraints on the transport properties of 
the system formed in high energy heavy-ion colli- 
sions. 

Experimentally it has not been possible to sepa- 
rate the contribution from elliptic flow (^2), fluctu- 
ations and non-flow from the data 0- This is due 
to the lack of knowledge of the probability distribu- 
tions for flow fluctuations. The fluctuations in V2 
can arise due to fluctuations in the eccentricity of 
the overlap region of the two colliding nuclei or de- 
viations of the participant plane from the reaction 



plane. While non-flow effects are those correlations 
among particles that are not related to the reaction 
plane (e.g due to resonance decay and jets). 



Asymmetric heavy-ion collisions could provide 
density profiles that are different or not accessible 
through symmetric heavy-ion collisions. For exam- 
ple a mid-central Cu on Au collision could lead to 
the Cu nucleus being occluded in the Au, leading to 
non zero odd harmonic of flow and being highly sen- 
sitive to early time dynamics. In addition asymmet- 
ric heavy-ion collisions could provide opportunity to 
disentangle the path length and energy dependence 
of parton energy loss. While in symmetric systems 
both initial energy density and transverse size of 
the medium increases with collision centrality, for 
asymmetric heavy-ion collisions it is possible to en- 
counter situations where the transverse dimension of 
the medium is same but one can get variation in en- 
ergy density. Keeping these aspects in mind, RHIC 
has now proposed to carry out asymmetric heavy-ion 
collisions program |8J] . In this report we only concen- 
trate on the study of the variation of 62, £3, their fluc- 
tuations, V2 and V3 for Au-I-Au, Au-I-Ag, Au-I-Cu, 
Cu-|-Cu and Au-|-Si collisions at ^/suN = 200 GeV 
using a multi-phase transport model (AMPT) with 
default settings. The goal being to simulate the ex- 
pectation of above observables in asymmetric heavy- 
ion collisions relative to symmetric heavy-ion colli- 
sions. 



The paper is organized as follows. The next sec- 
tion deals with the deflnition of €2, £3, their fluctu- 
ations, V2 and W3, along with a brief description of 
AMPT model. Section HI presents the results from 
the model calculation for the observables mentioned 
above. Finally we summarize our flndings in section 
IV. 
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II. AMPT MODEL AND DEFINITIONS 

The current work is based on the AMPT model 
with default settings |9!j . It uses the same initial con- 
ditions as considered in the Heavy Ion Jet Interac- 
tion Generator (HIJING) event generator HI- 
JING is a perturbative QCD inspired model which 
produces multiple minijet partons, these later get 
transformed into string configurations and then frag- 
ment to hadrons based on the Lund jet fragmenta- 
tion model fTT|. A Glauber model prescription is 
followed for obtaining the number of participating 
nucleons (A^part)- In AMPT, the minijet partons are 
made to undergo scattering before they are allowed 
to fragment into hadrons. These interactions could 
give rise to anisotropy in particle production along 
azimuthal direction. The event plane in AMPT is 
along the x-axis. There exists another version of 
the model, not used in this paper, called the string 
melting version where partonic interactions are con- 
sidered p]. 

We have followed the notations for the various ob- 
servables as studied in Ref [1]. The participant ec- 
centricity is defined as: 
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FIG. 1: (Color online) (a) Average eccentricity ((£2)) 
and (b) average triangularity ((£3)) as a function of 
Apart) for the various heavy- ion collisions at y/swN = 
200 GeV from AMPT model. 



(1) 



where r and (?!)part are the polar coordinate positions 
of participating nucleons in the AMPT model. ■02 is 
the angle of the minor axis of the ellipse defined by 
this region and is given as 



"02 



atan2 ((r^ sin(20part)) , {r'^ cos(20pa,.t))) + ti" 



(2) 

V2 which is the 2nd Fourier coefficient of the particle 
distribution with respect to '!/'2 and is given by 



V2 = (cos(2(0-^/>2))). 



(3) 



Similar to the definition of the eccentricity and 
elliptic flow, the participant triangularity, £3, and 
triangular flow, W3 are defined as: 



£3 



\J (r^ cos(3(/)part))^ + (r^ Sin(30part))' 



V3. = (cos(3((/) - Vs))) 



(4) 
(5) 



where ■03 is the angle of the minor axis of participant 
triangularity and is given by 



-03 



atan2 ((r^ sin(30part)) , {r'^ cos(30part))) + 



(6) 



III. RESULTS 

Figure [T] shows the (^2) and (£3) for asym- 
metric heavy ion collisions (Ag-|-Au, Cu-|-Au and 
Si+Au) compared to symmetric heavy-ion collisions 
(Au+Au and Cu+Cu) at ^siS" = 200 GeV as a 
function of A^part- We find for a given A^part, the (£2) 
values are higher for larger colliding ion sizes. How- 
ever for the low mass number colliding ions (Cu+Cu 
collisions and smaller) they are similar. No such 
large differences are observed for (£3). This suggests 
that within the framework of the AMPT model, the 
asymmetric heavy ion collisions can be used to con- 
strain the models dealing with second harmonic flow 
coefficient and its fluctuations. However, such col- 
lisions may not be that sensitive to studies dealing 
directly with triangularity and triangular flow. 

Figure [5] shows the fluctuations in £2 and £3 
expressed as the ratio of the corresponding root 
mean square (rms) values to their average values 
for both asymmetric heavy-ion collisions and sym- 
metric heavy-ion collisions at -^snn = 200 GeV as a 
function of Npart- We observe that for mid-central 
collisions (60 < Apart < 250) the fluctuations in £2 
increases as the system size decreases and then sat- 
urates for colliding systems of size comparable or 
smaller to Cu+Cu. In contrast the variation in flue- 
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FIG. 2: (Color online) (a) Ratio of root mean square 
(rms) value of £2 to {£2) and (b) rms of £3 to {£3) for 
various heavy ion collisions at y^SNN = 200 GeV using 
AMPT model. 
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FIG. 3: (Color online) {£^)/{£^)^, with n = 2,3, versus 
A'part for various heavy-ion collisions at ^snn = 200 GeV 
using the AMPT model. 



tuation of £3 as a function of system size is very 
small. 

Similar conclusions are obtained using a slightly 
different observable as proposed in Ref. pjj. The 
results of which from the AMPT model are shown 
in Fig. 131 It has been shown that the relative mag- 
nitude of Wn{4} and Wn{2} depends on the event-by- 
event fluctuations of if the non flow effects are 
small. Where the {4} and {2} represents the four 
particle and two particle cumulant methods to ex- 
tract flow coefficients respectively. It has been shown 
in Ref. ^ that 



Vn{2} 



e„{2} " 



(7) 



Figure m shows the same results as in Fig. [3] but as 
a function of fraction of collision centrality. Where 
corresponds to the most central collisions and 1 
corresponds to the most peripheral collisions. In ad- 
dition we also show for comparison results from a 
Glauber model simulation for Pb-|-Pb collisions at 
=2.76 TeV from Ref. [ij. Since the collid- 
ing ion species have different maximum A'part values, 
it is desirable to study how the fluctuations in £2 
and £3 vary for the case of same fraction of collision 
centrality. From Fig. |4] one observes that although 
the fluctuations in £3 is still similar for all the col- 
liding ion species studied, those for £2 now starts 
showing the ion size dependence more clearly. For a 
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FIG. 4: (Color online) (£^)/(£n)^i with n = 2, 3, versus 
centrality for various heavy-ion collisions at ^/swn = 200 
GeV. The Pb-I-Pb results corresponds to Glauber model 
simulations from Ref. fl^ at y^SNN = 2.76 TeV. 
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FIG. 5: (Color online) (^2) (a) and (113) (b) of charged 
particles at midrapidity for 80 < A'^part < 120 versus (£2) 
and (ea) respectively. The results shown are for various 
species of heavy-ion collisions at y^sjm = 200 GeV from 
AMPT model. The errors shown are statistical. 



FIG. 6: (Color online) V2 (a) and V3 (b) for charged 
particles as a function of transverse momentum (pt) at 
midrapidity for 80 < A'^part < 120 for various colliding 
systems at y/sjm ~ 200 from AMPT model. The errors 
are statistical. 



given fraction of collision centrality the fluctuations 
in S2 starts to increase as the colliding ion species 
size decreases, except for the most central and pe- 
ripheral collisions. These results show that a sys- 
tematic study of heavy-ion collisions with colliding 
system sizes ranging between Cu-|-Cu and Au-|-Au 
may provide a better handle on understanding £2 
fluctuations and its relevance to V2 fluctuations in 
the experimental measurements. 

Collective flow is mostly driven by the initial spa- 
tial anisotropy, hence it is expected that V2 and V3 
should be proportional to £2 and £3 respectively in 
high energy heavy-ion collisions. Figure [5] shows (W2) 
versus (£2) and {vs) versus (£3) from the AMPT 
model at midrapidity for 80 < A'part < 120 in heavy- 
ion collisions studied at ^/suN = 200 GeV. We ob- 
serve that indeed (W2) and (vs) are proportional to 
(£2) and (£3) respectively for all the collision sys- 
tems studied. However for a given value of (£2), 
the observed (W2) decreases with a decrease in the 
colliding system size. This could be due to larger 
fluctuations in £2 for the smaller systems compared 
to the larger colliding systems. However there is no 
such observed differences for (v^) versus (£3). This 
is consistent with not much difference in fluctuations 
of £3 for various colliding ions. 

Finally in Fig. EJa) we show the V2 for charged 
particles at midrapidity for mid-central collisions as 



a function of transverse momentum (pt) for vari- 
ous colliding systems at ^/swn = 200 GeV from the 
AMPT model. For the centrality range studied we 
see a clear dependence of V2 on the colliding system 
size. It increases with an increase in the colliding ion 
size. This is consistent with the results on variation 
of (112) with (£2) shown in Fig. [S] and the fluctua- 
tions in £2 shown in Figs. |3] and 21 In Fig. [BJb) the 
V3 versus shows a much smaller dependence on 
colliding system size compared to the above case. 

IV. SUMMARY 

In view of the proposed asymmetric heavy-ion col- 
lision program at RHIC in 2012, we have presented 
an AMPT model based study of £2, V2, £3, W3 and 
fluctuations in £2 and £3 for various asymmetric col- 
hding ion species of Ag-|-Au, Cu-|-Au and Si-|-Au at 
y/sNN = 200 GeV. These results are presented as 
a function of the number of participating nucleons, 
fraction of collision centrality and transverse mo- 
mentum. They are compared to results from sym- 
metric colliding systems of Au-t-Au and Cu-|-Cu col- 
lisions. 

Measurable signals of each of the above observ- 
ables are found for all the colliding systems studied. 
We find that while £2 and its fiuctuations, for a given 
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number of participating nucleons or fraction of colli- 
sion centrality, are highly dependent on the colliding 
ion type those for £3 and its fluctuations are very 
similar for all the colliding species studied. These 
results are reflected in the experimental observables 
such as V2 and W3. For the same £2 at midrapid- 
ity and for mid-central collisions, the proportional- 
ity constant between V2 and £2 seems to depend on 
the colliding system size. On the other hand, V2 vs. 
Pt decreases as the colliding system size decreases 
for the collision centrality range studied. However 



no such large sensitivity to colliding ion type on 
U3, £3 and fluctuations in £3 are observed from our 
AMPT model based study. Our study thus indicates 
that asymmetric heavy-ion collisions can be used to 
constrain models dealing with flow fluctuations in 
heavy-ion collisions but with greater sensitivity for 
V2 related observables than for V3. 
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